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°GMAT Basics (D. Conway) 

¢ Mission Design Walk-Through (S. Hughes) 
°GMAT Navigation Example (D. Conway) 

¢ CSALT Demo (S. Hughes) 

¢ Wrap-Up (S. Hughes) 
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|. GMAT Overview 


IV. Tour of the Script 
ll. Key Concepts Language ‘i 


Two Parallel Interfaces 


a a. Basic Syntax 
b. Resources and Commands 
c. Fields and Parameters b. eee Seen 
d. Execution Model c. Using Mat 
d. Using Parameters 
Ill. Tour of the Graphical © S°ver | 
User Interface eds 
Aeciiicontrale g. Best Practices 
b. Resources Tree : 
c. Mission Tree V. A Quick Demo 
d. Output Tree 
e. OrbitView 
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GMAT 
OVERVIEW 


¢ Design Started in 2002 
¢ Development Began in early 2003 
¢ Initial Fully Functional Build: 2004 
¢ First Public Demonstration/Release: 2007 
— Alpha/pre-Beta System 
— Used for Mission Design and Analysis 
¢ First Public Beta Release: R2011a 
¢ First Production Release: R2013a 
— First Operationally Ready Release 
— Target: The Advanced Composition Explorer (ACE) 
¢ Current Release: R2016a 
— First Production Ready Navigation Release 
— Target: The Solar and Heliospheric Observatory (SOHO) 
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System Features 


¢ Platforms Supported: Windows, Mac, Linux 
¢ External Interfaces: MATLAB and Python 
¢ Development Approach: Modified Open 
source 
— Developed Behind a NASA Firewall 
— Periodic Public Releases of Builds and Code 
— Supports a Robust Plugin Framework 
¢ Extensively Tested 
— More that 13000 Core Code Tests Run Nightly 
— More that 3000 GUI Tests 
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¢ Orbit design, optimization, 
and selection 

¢ Control design 

¢ Visualization 


¢ Orbit product generation 
and delivery 


¢ Event detection/prediction 


¢ Fuel bookkeeping & 
lifetime analysis 


¢ Propulsion system sizing 


GMAT Software Demonstration 


¢ Launch window analysis 
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Sensitivity and Monte Carlo 
analysis 


Navigation data simulation 
Orbit determination 


Maneuver planning and 
calibration 


Maneuver Support and 
reconstruction 


End-of-Life modelling 
Ephemeris prediction 


World-class quality software 


— TRLQ, Class B, (Part of Center-wide CMMI 
Accreditation) 


— Over 16,000+ automated script and GUI tests 


Large system with extensible design 
— 540k C++ LOC Core 
— Script, GUI, and plugin interfaces 


— 2 Interfaces to external systems (MATLAB and 


Python (under development) 


— 890k LOC from other libraries (GNOPT (Stanford 
Business Software). SPICE (JPL NAIF), Wx- 


Widgets, VF 13ad (Harwell), TSPlot Plotting 


Package (Thinking Systems, Inc.), Mars-GRAM 


model (MSFC) 
Enterprise level support 


— Large online support site (wiki, forums, issue 


tracker, downloads, etc) 


— Extensive Documentation (~1000 page User Guide 


and Reference Manual and ~100 pages of step-by- 


step tutorials) 


— Training (full-day live training courses and recorded 


training available via YouTube channel) 
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® General Mission Analysis Tool (GMAT) 


etx 


Window Help 


®@ DefauttOrbitview 


utput 
@ Defaultorbitview 


WS) ce S| > WH Be) 4 | © @ | animationéy © OJ G| curscript sync status: EAMONN 
DT: 
ion 


es of=)s fsa) | 


{@@ DefaultcroundTrackPlot 


B. Earthm2000Eq 

dA. Earthy2000Ec 

4. EarthFixed 

(8. EarthiCRF 
© Functi 


GMAT Wiki Home 


9 NASA missions 
¢ 5+ Discovery proposal efforts 


¢ 15 domestic and international 
universities 


° 6 OGAS 


¢ 12 contributing commercial firms 


A, @) 


THE UNIVERSITY 


| OF ARIZONA. a.i. solutions 


¢ 13 commercial firms using in open 


literature j rs 
¢ 30+ independent peer reviewed Universitat oY OLEAN = 
publications citing analysis hia al JPL 
performed using GMAT Aly «= eC S a sordnligeor opel oa neem 
GMAT is used y A]|M WS 


European Space Agency 


world-wide 
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fp | 
Ames Research Center 


AAS Guidance and Control Conf 
Feb. 7, 2017 
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KEY 
CONCEPTS 


OrbitView Window Ground Track Window Graphics Window 


® GMAT - General Mission Analysis Tool 


Menu Bar — 
| File Edit Window Help 
Toolbar | ey EY | f | PH BS! SC! OC O | animdtton BG DB G | curscipt sync status: 
R 
Resources Tab | ee © DefaultOrbitView 
=|) Spacecra f " . an 
— DefaultSC | ? Coma Beenicés 
Mission Tab © Formations mit 
{} Ground Station 
{© Hardware 
Output Tab a 
6%, DefaultiB 
=} Propagators 
% DefaultProp 
Resources Tree 4 SolarSystem 


+) Solvers 
=) Output 
@ DefaultOrbitView 
tad DefaultGroundTrackPlot 
=) Interfaces 
{m Matlab 
{5 Scripts 
{© Variables/Arrays/Strings 
=) Coordinate Systems 
di. Earthm32000Eq 
i. Earthm32000Ec 
i. EarthFixed 
ct. EarthICRF 
{ Functions 


Message Window Running mission. .. 


Mission run completed. 
===> Total Run Time: 0.661000 seconds 


status RP Gonware Demonstration 


°GMAT is like MATLAB: 


— You write a program (a “mission’), then run it to 
generate output 


¢ Not like Excel 


— Cannot generate output or manipulate results 
without rerunning 
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¢ Batch execution model 
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® C:\Users\jjparker\Desktop\MySd% 
@ File Edit Window Help 


| iar BI is oo e\>r ug 
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Script 


General Mission Analysis Tool(GMAT) Script 
$Created: 2013-01-23 09:47:07 


Create Spacecraft DefaultSc; 

GMAT DefaultSC.DateFormat = TAIModJulian; 

GMAT DefaultSC.Epoch = '21545'; 

GMAT DefaultSC.CoordinateSystem = EarthMJ2000Eq; 
GMAT DefaultSC.DisplayStateType = Cartesian; 
GMAT DefaultSC.X = 7100; 

GMAT DefaultSC.Y = 0; 

GMAT DefaultSC.Z = 1300; 

GMAT DefaultSC.VX = 0; 
GMAT DefaultSC.VY = 7.35; 
GMAT DefaultSC.VZ = i; 

GMAT DefaultSC.DryMass = 850; 
GMAT DefaultSC.Cd = 2.2; 


GUI and script are nearly interchangeable (but not totally). 


GMAT Software Demonstration 
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Resources Commands 


¢ Participants in a GMAT ¢ Events ina GMAT mission 
mission - Represent the actions 

¢ Represent the “things” that taken on the resources 
will be manipulated * Think of them as methods 

¢ Think of them as objects, or functions 


with properties 
¢ Most are “fixed” when the 
mission starts 
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Fields 


¢ Properties you can set on 
a resource 

¢ Examples: 
— Spacecraft.Epoch 
— Thruster.DecrementMass 
— ReportFile.Filename 


Parameters 


° Properties you can 
calculate during the 
mission 

¢ Parameters often have 
dependencies 


¢ Examples: 
— Spacecraft.Earth.Altitude 


— Spacecraft. EarthMjJ2Q00Eq 
.BVectorAngle 


¢ Sometimes a property is both a field and a parameter. 
¢ Examples: Spacecraft.SMA, FuelTank.FuelMass 
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GMAT Software Demonstration 


TOUR OF THE 
GRAPHICAL 
USER 
INTERFACE 


GMAT Software Demonstration 


¢ Contains all configured 
resources in the mission 
¢ Grouped into folders by type: 
— Spacecratt 
— Hardware 

— Burns 

— Output 

— SolarSystem 
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¢ Contains the Mission 
sequence—sequence of all 
configured commands 

¢ Special features: 
— Docking & undocking 
— Filtering controls 
— Command Summary 


GMAT Software Demonstration 


SWF BD | wR © Be | ay << <e |(F] 
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¢ Contains all outout products 


¢ Populated after mission 
execution 


GMAT Software Demonstration 
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¢ 8D graphics window 
¢ Most complex of the graphical output types 
— Others include: XYPlot (2D plotting), Ground TrackPlot 
(2D mapping) 
¢ Mouse controls: 
— Left button: rotation 
— Right button: zoom (horizontal motion) 
— Middle button: rotation normal to screen 
¢ Configuration includes: 
— Camera controls 
— Resources to draw 
— Visual elements 
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TOUR OF THE 
SCRIPT 
MAN eley\e) = 


* Syntax is based on MATLAB 

¢ Single-line statements w/ optional line 
continuations 

¢ Case sensitive 

¢ Loosely typed 

¢ Begin/End block statements 


¢ Resources are created before used (except 
special defaults like SolarSystem) 


AAS Guidance and Control Conference, DA 
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¢ Script is divided into two sections: 

— Initialization (at the top) 

— Mission Sequence (at the bottom) 

— Divided by the BeginMissionSequence command 
¢ Initialization -> Resources Tree 

— Static assignment only 
¢ Mission Sequence -> Mission Tree 


— Manipulation of existing resources, cannot create 
new ones 
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Create Spacecraft sat 
Ssat.SMA = 7000 


Create ReportFile r 
r.Filename 


"MyReport.txt' 


BeginMissionSequence 
Report ‘Write SMA' r sat.SMA 


GMAT Software Demonstration 
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Ee Solvers 


¢ Math syntax is based on MATLAB 
¢ Operators are matrix-aware 
| Built-in F 


- ig cos 
[es add asin 
- subtract aeee atan 
* multiply atan2 log 
/ divide se SEP 
' transpose DegToRad RadToDeg 
‘power abs sqrt 
norm det 
inv | 
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Create Spacecraft SC 

SC.SMA = 7100 

Create Variable period, mu, pl 
mu = 398600.4415 


BeginMissionSequence 


pi = acos(-1) 
period = 2 * pi * sqrt(SC.SMA%3/mu) 
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¢ Parameters can have one of two types of 
dependencies (or neither): 
— Central body 
— Coordinate system 

¢ They are calculated on the fly when they 
are used: 
— Spacecraft.MarsFixed.X 
— Spacecraft.Earth.BetaAngle 


- If omitted, default dependency is used 


ee : AAS Guidance and Control Conference, 
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Cr ets U f 


Create Spacecraft SC 
SC.CoordinateSystem = MarsFixed 
Create ReportFile r 
BeginMissionSequence 


% using parameters 
Report r SC.EarthMJ2000Eq.X 
Report r SC.Earth.BetaAngle 
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¢ Three control flow statements: 
— |f/Else — execute If a conditional is true 
— While — loop while a condition is true 
— For — loop a certain number of times 


If SC.Earth.Altitude < 300 
% dO a maneuver 
Else 


% continue 
EndIf 
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¢ Three types of solvers: 
— Targeter (using Differential Corrector) 
— Optimizers (using either optimizer) 
— Estimator (using Batch Estimator and, soon, EKF) 
¢ Similar to loops, with specific nested 
commands: 
— Target: Vary, Achieve 
— Optimize: Vary, NonlinearConstraint, 
Minimize 
¢ See the tutorials for examples 
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¢ Performing a Hohmann Transfer 
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DISCUSSION 
AND 
QUESTIONS 
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GMAT Application to GSFC Mission Design 
Steven P. Hughes 
08 Feb. 2016 


This presentation is based on presentations provided by the TESS project 
and used with their permission. Author attributions are listed throughout. 


General Mi 
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Sess : Agenda 


me a a 


Mission Overview 
Requirements 

Trajectory Design Process 
Other Areas 


Solution Generation Process 


Finite Burn Modeling We won’t have time to 


Launch Vehicle Dispersion Analysis discuss these but all of 
these activities were 


Maneuver Planning 
performed in GMAT 


Launch Window Analysis 


Se lee ee 


Conclusions 


NOTE: This is a snapshot from TESS PDR and some details have changed. 
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&®ss Transiting Exoplanet Survey Satellite 


a 


Bama 01: Mission Overview 
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Primary Goal: Discover Transiting Earths and Super- 
Earths Orbiting Bright, Nearby Stars 


= Rocky Planets & Water Worlds 
= Habitable Planets 


Star + Planet Dayside 
pe 


Discover the “Best” “1000 Small Exoplanets 


- = “Best” Means “Readily Characterizable” 
Star Planet Nightside . Bright Host Stars 
7 Transit 


Z * Measurable Mass & Atmospheric Properties 
Star — Planet Shadow 


= Present: Only 3 small transiting exoplanets orbiting 
bright hosts are known 


Large Area Survey of Bright Stars 
= F, G K dwarfs: +4 to +12 magnitude 
= M dwarfs known within “60 parsecs 


= “All sky” observations in 2 years: 
e > 200,000 target stars at <2 min cadence 
¢ > 20,000,000 stars in full frames at 30 min cadence 
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_SWST ome, 
continuous %% 
viewing zone *** 


Anti-Solar segments Coverage of ecliptic poles drives 
drive +/- 15 deg Pitch angle (nominally 54 deg) 


¢ Concentration of coverage at the ecliptic poles for JWST. 


¢ Sacrifice of coverage in the ecliptic because Kepler-2 is already mapping that 
region. 
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Launch to Science Orbit Timeline 


© Burn if necessary 
@ DvVBurn 


ToM @  Perigee Passage 


Phasing loops 2 & 3 
variable by launch date r, = 376,300 km (59 RE) 


Lunar 
Swing-by 


ra=250,000km At 


Cal 


s 


p= 108,400 km 
(17 RE) TM 


¢----------------->, 


a, = 200 km 
Injection 
Transfer Orbit Sey Science Orbit Science Orbit 
1 2 
Ascent and Commissioning (60 days) Science Operations 
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{ess Nominal Aug 10 solution: Inertial frame 
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Nominal Aug 10 solution: Rotating frame 


For a loop in 
the1st quadrant, 
the Moon is 
behind and lowers 
perigee 


For a loop in the 
4'h quadrant, the 
Moon is ahead 
and raises 
perigee 
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Roughly 47 degree plane change at flyby 


Earth Inertial Axes 


5 Oct 2017 01:19:00.000 Time Step: 1800.00 sec 
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Sess 


Key L2 Mission Design Requirements 


1B) Title Requirement Summary 

MRD_2 Mission Life 2-year mission + 2-month commissioning 

MRD_10 Observation Period HASO duration = 12.5 days per orbit 

MRD_54 Launch Period Launch opportunities on at least 5 days days per lunar cycle 

MRD_55 Launch Window 30-Second Launch window 

MRD_42 Ascent and Commissioning Duration | Achieve mission orbit within 2 months after launch 

MRD_51 Mission Orbit 2:1 lunar-resonant orbit 

MRD_52 Maximum Range in LAHO Perigee < 22 Re 

MRD_101 | Mission Maximum Range Apogee < 90 Re 

MRD_53 Avoidance of Geosynchronous Orbit | Orbit does not intersect GEO band for mission + 100 years 
(TBD) 

MRD_56 Eclipse Frequency and Duration No eclipses longer than 5 hours and not to exceed 14 in 
number (duration = umbra + 0.5*penumbra 

MRD_104 | Delta-V Allocation Total AV < 215 m/s (99% probability) 

MRD_129 | Longest Single Maneuver Longest continuous maneuver < 95 m/s 

MRD_85 Sun in Instrument Boresight FOV exclusion of 54°x126° (TBR) for 15 minutes (TBR) 

MRD_64 Vista Nchi cae arava see eee 
rraneuver(Deleted} 


Green Requirements are Orbit Design Drivers 
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Consistent with EXP-TESS-GSFC-RQMT-0001 Rev B 
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&®ss Transiting Exoplanet Survey Satellite 


| i 03: Trajectory Design 
One Process 
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Sess Overview 


The TESS trajectory design process is based on three components: 


« Theoretical basis 
= Kozai constant 


= Tisserand condition 


¢ Two-body patched-conic first guess 


= Implementation of theory to approximate final trajectory 


¢ High-fidelity targeting 


= Transitions approximate first guess to realistic final solution 
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apss Implementation Overview 
¢ General Mission Analysis Tool (GMAT) used for t= 
implementation of design ||. 
= GSFC’s in-house high-fidelity trajectory design software Sa 
¢ Uses first guess to seed numerical targeting Vea 
algorithm 


approx. 
flyby 
epoch Two-Stage 
High-Fidelity erties 
Targeting 


approx. 

phasing 
loop 

duration 
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&bss GMAT Design Approach 


¢ Two targeting stages 


¢ Stage 1: Design from Translunar 
Injection (TLI) through flyby to 
Science Orbit 
= Multiple-shooting process 
= Starts with patched-conic first guess 


¢ Stage 2: Backwards design from 
converged mission orbit to launch 
vehicle separation (adding 
phasing loops) 
= Single-shooting process 


= Starts with converged outbound 
solution + 2-body phasing loops 
guess 
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Sess Outbound Sequence Overview 


¢ Multiple-shooting approach w/ 5 segments 


® ( 


TLI ‘ : swingby = : PAM , 


PP4 
Xto\" = X 


swingby tswingby Ssmtaalny = Xpam 
swingby 


@ control point 
O patch point 


¢ Start with patched-conic initial guess for each segment 


¢ GMAT targeting sequence used to find smooth solution from 
segmented initial guess 
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ass . Implementation Overview 


¢ General Mission Analysis Tool (GMAT) used for 


implementation of design 
= GSFC’s in-house high-fidelity trajectory design software 


¢ Uses first guess to seed numerical targeting 
algorithm 


approx. 
flyby 
epoch Two-Stage 
High-Fidelity erties 
Targeting 


approx. 

phasing 
loop 

duration 
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¢ The TESS trajectory has two critical features: 
= Transfer orbit (result of lunar flyby) 


= 2:1 lunar resonant mission orbit 


TESS HEO Orbit Definitions 


Lunar Resonance Phasing (LRP) or PLEP Angle Transfer Orbit 


“Leh @ Lunar Fy 


@ TESS PAM Burn 


cosAOP =a _ar(l-er') dy 
er r 


| Earth-Moon Plane; Descending Outbound 
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&bss Tisserand Criterion 


¢ The Tisserand criterion holds that a quantity T is constant before and after 
a flyby: 
1 
= Fn + cos(i) Ja(1 — e?) 


« Here ais semimajor axis (scaled by distance between the primary bodies), 
e is eccentricity and tis inclination to the orbit plane of the primaries 


¢ The Tisserand criterion is used for TESS to design the lunar flyby. 


= We choose the value of T to obtain the desired orbit properties of the 
transfer orbit after flyby to mission orbit. 

= The transfer orbit shape is driven by a timing condition: the need for the 
spacecraft at Post Lunar Encounter Perigee (PLEP) to nearly line up with 
the Moon. The spacecraft-Earth-Moon angle at perigee is called PLEP 
misalignment or the Lunar Resonant Phase Angle. 


= We then use the value of T to infer the shape of the orbit before flyby 
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Sess Kozai Mechanism 


¢ The Kozai Mechanism describes the long-term evolution of a 
highly eccentric, highly inclined orbit due to a third body 
(Moon). 

¢ The Kozai model implies that: 


= Orbit semimajor axis is conserved 


= Kozai parameter K = cos(i)vV1 — e? is constant, where e is eccentricity and i is 
inclination to the Moon orbit plane 


¢ Kozai mechanism predicts 


= Eccentricity and inclination oscillate in unison, with a period of about 8 
years for a TESS-like orbit. (Therefore, perigee radius and inclination 
oscillate together.) 


= AOP relative to the Moon librates around 90 deg or 270 deg, if the initial 
inclination is higher than critical inclination 39.2 deg 
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Sasetine-Tesal 7 - 17 Jum 2016 17:99:12 


orbit period (day) 


— 2020 2025 =e 2030 2035 


Orbit period oscillates around ideal 13.65 days with amplitude 
near 1 day 
Thus we do not need to start at ideal resonance to achieve 
resonance on average 

¢ We change mission orbit period by changing the Period Adjust 


Maneuver (PAM) at mission orbit insertion 
IBEX extended mission design used this feature to select an initial 


orbit that achievedJquaiechnes ducakispas tar saveiat years 
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Sess Kozai Mechanism (cont’d) 


¢ Kozai mechanism is relevant to TESS because 
= We want mission perigee radius to remain between 6.6 Re (GEO) and 22 Re 
= We want mission ecliptic AOP to remain near 90 deg or 270 deg, so line of 
apsides stays out of ecliptic plane, and so long eclipses cannot occur near apogee 
¢ For TESS orbit, e = 0.55so0 K = 0.65 implies i = 39 deg 


¢ We exploit the fact that the lunar plane and ecliptic plane are near the 
same, only 5 deg apart. 

¢ Perturbing forces (especially the Sun) imply that the Kozai mechanism does 
not work exactly in the full force model. Nevertheless, like CR3BP, the Kozai 
mechanism is a useful technique for orbit design 


Methods described by Aerospace Corp in CSR and flight dynamics paper “A High Earth, Lunar Resonant 
Orbit For Lower Cost Space Science Missions” by Gangestad, Henning, Persinger and Ricker (AAS 13- 
810) 
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caauiiciaiaiae : Evolution of perigee 
radius (green) and lunar 


«2 5 5 . 
2” bed inclination (red) over 20 
: = years. The oscillation 
2 += 20 me . . 
2 36. {3 period is about 10 years 
2 a6 & 
2 30 + 2 


20200 = (2025 ms 2030 = (2038 


Savesine-Tesst T . 17 Jun 2016 18:00:51 


‘aiid 
Evolution of ecliptic AOP 
oe 4 over 20 years. 
> 90~: 
S 
< 
8 60 «= 


— 2020 2025 saliniad 2030 2035 


See also: Dichmann, Parker, Williams, Mendelsohn: Trajectory Design for the Transiting 
Exoplanet Survey Satellite. ISSFD 2014 
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apes Implementation Overview 
¢ General Mission Analysis Tool (GMAT) used for t= 
implementation of design ||. 
= GSFC’s in-house high-fidelity trajectory design software Sa 
¢ Uses first guess to seed numerical targeting Vea 
algorithm 


approx. 1 
flyb 
oh Two-Stage 

High- End-to-End 
approx Fidelity Solution 

; r ti 
phasing al ia 
loop 

duration 
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Sess GMAT Design Approach 


Two targeting sequences 
Stage 1: Design from Translunar 
Injection (TLI) through flyby to 
Science Orbit 
= Multiple-shooting process 
¢ Stage 2: Backwards design from 
converged mission orbit to launch 
vehicle separation (adding phasing 
loops) 
= Single-shooting process 
= Starts with converged outbound 
solution + 2-body phasing loops guess 
« Both stages use VF13 NLP solver as 
robust targeter 
= Seeks feasible solution only; not 
optimizing 
¢ Final 3™ stage: forward-propagation 
from SEP to check constraints 


Change since PDR Peer Review 
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ess Modeling Assumptions 


¢ All analyses share common force models, spacecraft 
parameters, solar system models, to the extent practical. 


Spacecraft model 


201.9 kg *Low dry mass estimate, used 
to model worst-case SRP 
: effect & kept for continuity. 
2 Current mass estimate is used 
in finite burn analysis. 


Force modeling Phasing loops Flyby Mission orbit 


JGM-2 40x40 | Moon JGM-2 8x8 
point mass 


Solar system ephem 


DE421 
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ass Stage 1: Outbound Sequence Constraints 


Value Description 


Fixes TLI at approximate LV insertion inclination 
600 km Phasing loop perigee altitude 
oa Fixes TLI at perigee 


Mission orbit perigee 17 RE Design value for min/max perigee behavior 
radius 


Fixes PAM at perigee 


Maximum misalignment from resonant condition 


Mission orbit energy 2:1 resonance Energy from SMA consistent with 2:1 resonant 
condition 
Mission orbit Kozai 0.60<K<0.80 | Controls long-term perigee behavior 


Position/velocity Position/velocity continuity between all segments 
continuity 
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Outbound Sequence Overview 


Guess with Discontinuities Solution, No Discontinuities 
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Phasing Loops Sequence Overview 


Starts with converged outbound solution 
Back-propagates from PAM through flyby to TLI 
Uses targeting sequence to add on phasing loops 


Two-body initial guess for A1—A3, P1—P3 burns 


Insertion constraint is now enforced at insertion, not at TLI 


Small out-of-plane components are added to PAM to correct inclination 
at TLI 


This is a side effect of the two-stage approach; would go away in an 
end-to-end solution 
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ass Stage 2: Phasing Loops Constraints 


Parameter Value DT =XYol gi eldcey a 


P1—P3 altitude >600 km Phasing loop perigee altitude 


Separation altitude 200 km LV requirement 
Separation inclination 28.5° TOD LV requirement 


Separation epoch match launch modeling & Analytical model based on launch site 
desired phasing loop duration 


275,000 km RE 
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Final Converged Solution 
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GMAT Navigation 


Darrel J. Conway, Thinking Systems, Inc. 


This presentation Uses Materials in the GMAT Simulation and 
Estimation Tutorials, Written by D. S. Cooley of NASA GSFC 
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l. GMAT Navigation Overview 
ll. Estimation Walkthrough 
lll. Discussion and Questions 
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GMAT 
NAVIGATION 
OVERVIEW 


¢ Estimators 

— Batch Least Squares 

— Under Development: Extended Kalman Filter 
¢ Measurement Simulator 
¢ Measurement Models 

— DSN 2-way Range and Doppler 

— (Alpha) GN 2-way Range and Doppler 

— (Alpha) SN 4-Leg Range and 5-Leg Doppler 
¢ Solve-For Parameters: 

— Cartesian State 

— CdandCr 

— Measurement Bias 
¢ Built as a GMAT Plug-In 

— Navigation Features Script Based 

— Alpha Measurements Disabled by Default 
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Bet epee ae ee : AAS Guidance and Control Conference, 

GMAT Software Demonstration : Feb 7 O17 aoa 68 
rT ~ f/f, CULE 


¢ Coded from Tapley, Schutz and Born, 
Statistical Orbit Determination 
¢ Extended with Features from GIDS 
— Robust Inversion Code 
— Sigma Editing 
— Detailed Estimation Output 
— GTDS Based Convergence Criteria 
¢ Provides Results Comparable to GIDS 
¢ Approved for Operations for SOHO 
— Pending successful ORR, 2017-01-31 


GMAT Software Demonstration 


¢ Measurements Contained in 
“TrackingFileSet’s 

- Specified by the Path 
Followed by Signals 

- R2016a Measurement Models: 
— DSN Range, in Range Units 
— DSN Doppler 

¢ Aloha Measurements: 
— 2-way Range and Doppler 
— 4-Leg Range and 5-Leg 

Doppler (TDRSS) 
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GMAT Software Demonstration 


4 
Solve-For 

qi 

a 


; 
| 

I 
: 


Hardware 


_ Media Corr’ns 


tg ce) mm \Vi(ele(s) 


q 


4 


~” 


GMAT Nav Components 


@ a » 4 » 
aay] Tracking Files -Conv'ce Criteria 
q . 4 
a — @ > | Sigma Editor — 
ce ee 
q 4 : 
» » 
q » d 
ce <? < 
S dj 
‘ 5 SolveFor | 
. © Bias, Sigma | 
S y % d 


ESTIMATION 
WALK 
THROUGH 


Antenna 


 Harcnare 


™N 
} 
| 


antenna 


Media Corr’ns | 


a 


_ Error Model 


ww 
“a 


~ Tracking Files 


| 
\ 


(/ 


| Solve-For 


| Taltcre[e-icele 
me) Cem (ele(=) | 


NX 


= » 
~ 


| Bias, Sigma | 


_Conv’ce Criteria — 


/ Sigma Editor | 
“Tracking File Set | 
-Propagator 


Va 


Bice File Set | 


| Propagator 


\ 


ff 
f 


SS 


| Solve-For | 


Hardware 


4 


~ 


Error Model Error Model 5 


| Media Corr’ns 


Hardware 


PNalc\alare! 


Antenna 


| Integrator 
Force Model — 


[ Bias, Sigma 


[ Solve-For | 


a») 


4 


_ Conv’ce Criteria | 


_ Sigma Editor 
Tracking File Set | 
-Propagator 


\ y 
SS i 
4 a = = ~~ 


“Tracking File Set | 


- Propagator 


S 


Configuring Resources ( 


" 


Solve-For | 


\ 
Xx 
SS 


Media Corr’ns | 


Hardware 


[ . Error Model i 


4 


Integrator 
Force Model 


Solve-For 


Bias, Sigma 


| Tracking Files | | 


X 


| Conv’ce Criteria | 


Sigma Editor 


| Tracking File Set | 
Propagator | 


| Tracking File Set | 
| Propagator 


X 


Solve-For 


Hardware 


y 


X 


Ww 


mtelaelirclas 


Media Corr’ns 


| =f c@) am \V(ele(>)| 


y 


= 


PNalccialarel 


PNalisialarel 


PNalc\alarel 


| Tracking Files 


“a 


| Integrator | 


[ Force Model 


X 


| Solve-For 


| Conv’ce Criteria | 

Sigma Editor 

“Tracking File Set 
-Propagator 


“Tracking File Set 


| Bias, Sigma 


I [ Propagator | 


¢One Command: RunEstimator 


File Edit Window Help 


KCAL 


Resources Mission Output 


-& Spacecraft 
VG Sat 


© \gsfegmat\GmatDevWin\zpplication\samples\Tut_Orbit_Estimation_using_DSN_Range_and_Doppler_Data.script - General Mission Analysis Tool (GMAT) 


|B HE ts SF | © © | mrirai 


{© Formations 
= Ground Station 
§é CAN 
6 Gps 
§é Map 
{© Hardware 
& Burns 
2 Propagators 
% Prop 


& 
8 
© 


° 
° 

° 

rs 


° 
2 
S 

% 


4 SolarSystem 272545 27258 
4G Solvers 
{ Boundary Value Solvers 
© Optimizers 


272615 27265 27268.5 27272 272755 


(>| Tut_Orbit_Estimation_usin 
3 Variables/Arrays/Strings 
(x) PosError 
(x) VelError 
49 Coordinate Systems 
dA. EarthMu2000Eq 
. Earth2000Ec 


fel) 


Oo GUY Script Sync Sat TE 


(eis) 


-13.6 


272 


CaF? 
oJ 


0 


G00 S00 
ogo G00 
oo BRO o 


° 
ao 
[10 0 Gore 


° 


‘ODSNsimData_{GDS,Sat,GDS}_DSNRange Residuals| 


27251.69 


27255.28 27258.87 


(27262.46 


27266.05 27269.64 27273.23 


[ODSNsimData_{MAD,Sat,MAD}_DSNRange Residuals| 


A. EarthFixed 272545 27258 

df. EarthiCRF 

-¢. SunMJ2000Eq 
Event Locators 


27261.5 27265 27268.5 27272 272755 


{ele lis 


© Functions Ss 
0.074E T T ° listenin arm aa Residuals] 
° ° ° o 

€o%o0 ry 

0037-8 BP OO OOS Og oP 
oe So of » BB oO ® 8 g 

of PF BIB LEK Egy VPS Fo 0 F Po 4 

9% oo g8 eo oe ee 6 Oo, 8° PY 
-0.037 - BP 5 o@ 2 o 8 @ 8 2o 08 24 

Oe a0 ° 
-0.074- o----0. ae J 
oa ! e 1 \ \ ! 
27251.69 2725528 © 27258.87 2726246 © -27266.05 + 27269.64 2727323 
« > 


|e 


[ODSNsimData_{CAN,Sat,CAN}_Doppler Residuals} 
v o 


27251.69 


27255.28 (27258.87 


(27262.46 


27266.05 27269.64 27273.23 


[e|e I 


° 
B25 we 
@--2. 
Seve gs, 8 
° 9 S 9 © 
eB oO 
gabe sees 
3 a 6 
Bois 5? 3S 


3 Bi" 
: cae 
She shee se | 
See se oF sR 
BF-8-8 yo $8.8 a 
8 O-'9 = 


‘ODSNsimData_{CAN,Sat,CAN}_DSNRange Residuals| 


27251.69 


27255.28 27258.87 


(27262.46 


27266.05 27269.64 27273.23 


Estimation converged! 
[1 - RMSP/RMSBI = | 1- 1.10756 / 1.10756] = 1.39866e-007 is less than RelativeTol, 0.0001 


Final Estimated State: 
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¢ Output Data: 


— Estimation Status (Converged, or Failure 
Reason) 


— Final Estimated State 


— Final Covariances and Correlations 
— Detailed Estimation Report 
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Fig. 1: Schematic of linkages for multiple-phase optimal control problem. The example 
shown in the picture consists of five phases where the ends of phases 1, 2, and 3 are 
linked to the starts of phases 2, 3, and 4, respectively, while the end of phase 2 is linked 


Phases 2 and 3 Connected 
Phases 2 and 5 Connected 


Pianes dando Connecied Phases 3 and 4 Connected 
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Phase 1 


Trajectory 


Phase'5 


to the start of phase 5. 


Note, the figure above was taken from Patterson, M, and Rao, A, A MATLAB Software for Solving Multiple-Phase Optimal Control Problems Using hp—Adaptive 


Gaussian Quadrature Collocation Methods and Sparse Nonlinear Programming”. 
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time 


or Stories: Trajectory 
iD [Description ~~=~S~C~*~<“‘<‘<“<CSCSCSCSS | 


CO.TD-1 Basic Orbit Transfer: | need to design an optimal maneuver to transfer to 
a nearby orbit. 

CO.TD-2 Body to body transfer: | need to design a trajectory that does a direct 
transfer from one body to another, in my case, from Earth to Moon. 

CO.TD-3 Interplanetary Multiple flyby: | need to design a trajectory that leaves 
Earth, performs multiple gravitational flyoys, and arrives at a final target 
celestial body. | have a guess at the flyby times, order, and flyby 
conditions at each body. 

CO.TD-4 Spiral: | need to design a trajectory that spirals into and/or away from a 
celestial body requiring many (100s) of revolutions. 


Libration point maintenance: | need to use collocation to improve the 
robustness and quality of libration stationkeeping solutions. 


GMAT Software 
Demonstration 


Ip |Deseription 


CO.TD-6 Weak stability transfer: | need to use collocation to find trajectory 
transfers in a weak stability, multi-body environment. 

CO.TD-7 Formation: | need to use collocation to optimize maneuvers of multiple 
spacecraft. 

CO.TD-8 Rendezvous/Prox Ops: | need to use collocation to optimize rendezvous 
and docking operations and proximity operations. 


CO.TD-9 Probe Separation: | need to use collocation to determine optimal 
trajectories for probe release and spacecraft separation planning. 
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GMAT High Level Design 
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Problem Name Truth Source GMAT Prototype External Truth (gmat-truth)/truth 
Optimal Cost Function Optimal Cost Function 


Rayleigh (control path constraint) | sos | 44.72093885 44.7209362 5.925E-08 
Rayleigh (state and path constraint) | sos | 44.8044450 44.8044433 3.794E-08 


HyperSensitive 6.72412985 6.72412325 9.815E-07 
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Interplanetary Application with Comparison to EMTG 
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A significant part of certification of navigation functionality is using 
GMAT in the FDF operational environment. 
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SDO Operations 
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We are performing final testing for using GMAT as the core tool in the SDO MOC for 
maneuver planning and product generation. 
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GMAT API and Potential Future Efforts BY 


MATLAB 


° (ee 


GMAT 
Onboard 


Flight Dynamics Library 


python” 


MATLAB 


& 


Java 


> 


We are currently in the early phases of developing a GMAT API, and considering the possibility of 
deconstructing GMAT to allow components to be used onboard for 
real time mission planning and opportunistic science. 
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Finding Out More 


e Main information portal: gmatcentral.org 


e For New Users 
= Obtain GMAT 
# Training Videos 
# Training Material 
= Sample Missions (scripts distributed with application) 
= User Guide (also pdf distributed with application) 


e For New Developers 
# Obtain Code 
# Style Guide 


" Compilation Instructions 
= Design Docs (also pdf distributed with application) 
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Thanks! Questions? 
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Past Release Summary Sy 


e GMAT R2013a = Customizable orbit segment colors 
= First production (non-beta) release = Mars-GRAM 2005 atmosphere model 
=" Focused entirely on QA and = LHS parameter dependencies 

documentation = New solver algorithms 
= Very few new features—but many 
improved e GMAT R2015a 
= New support for ICRF coordinate = GMAT Functions 
evens = Python Interface 

e GMAT R2013b (internal) = Eclipse Location 
= First operationally-certified release = Ground station contact location 
= Focused on ACE mission = SNOPT Optimizer 

requirements # Space weather modelling 


= 3D models for celestial bodies 


e GMAT R2014a = Solver status window 


= Public release of all R2013b features 
= State representations 
= Attitude models 
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Ongoing Navigation Development S& 


= 2009 - 2011 =» 2014-2015 
= Began evaluation of = Measurement model re-design 
GMAT as a possible based on GEODYN principles 
navigation tool in 2009 = User interface re-design for 
= Worked with AFRL and usability based on FDF feedback 
IRAD funding to design - Testing against flight data 
and implement a | ; 
navigation subsystem = Improved batch estimator 
and demonstrate = New data types 
feasibility. = Measurement editing 
= Key Conclusion: GMAT =| dR 
could perform OD without ee epemne 
significant design = Improved bias modelling 
changes. = Improved inverse algorithms for 
=» 2012-2013 normal equations 
: Interplanetary models = New Solve-fors 
dynamics models = Low thrust navigation studies 
= DSN data types = Major testing effort in FDF 


GMAT was selected as the core tool for GSFC navigation and is 
J) AAS Guidance a/Preparingdor operational use,in fall of 2016 100 


Extensibility 


e GMAT’s modern architecture was designed for 
extensibility 


e Extensible System Interfaces 
= MATLAB 
= Python 
= API under development 
= Plugins 


e Multiple User Interfaces 
= Script 
= GUI 
= Command line 
= API under development 


External Processes 
\ 


= wel 


GMAT Engine 


e Extensible model subsystems > 
= Dynamics Models 
= Environment Models > 
= Estimators 


= Measurements 
= Propagators 


python &, 
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L3 Mission Design Requirements vasa 


Requirement Complian 
for) 
L3_FD_1 | MRD_10 | Mission Orbit SMA | The target mission orbit Semi-Major Axis Comply. 
(SMA) shall be 38 Re. Design 
MRD_51 constraint. 


TESS HEO Orbit Definitions 


Lunar Resonance Phasing (LRP) or PLEP Angle 


@ TESS PAM Burn 
ar(l -er'), 


ine of Nodes 


cosAOP = Ja - 
er 


I Earth-Moon Plane; Descending Outbound 


Consistent with EXP-TESS- 


FC-ROMT-0015 Rev (- 
erence, Feb. 7, 2017 Coreen Oe ea “08 


L3 Mission Design Requirements 


S 


Parent Requirement Complian 
1B) ce 
L3_FD_3 | MRD_53 | Mission Orbit FD shall target a mission orbit with a minimum | Comply. 
Minimum Perigee perigee that shall stay above GEO radius + Results 
200 km. shown to 
100 years. 
L3_FD_2 | MRD_52 | Mission Orbit FD shall target a mission orbit with a Comply. 
9 Maximum Perigee | maximum perigee that shall stay below 22 Re | All <20.5 
for the duration of the mission. Re 
L3_FD_3 | MRD_10 | Transfer Orbit FD shall target a lunar flyby that results ina Comply. 
0 1 Maximum Apogee __| transfer orbit with a maximum apogee less All <80 Re 
L3| FD_{29, 30, 38} TESS HEO Orbit Definitions 
replace old L3 FD 3 in Lunar Resonance Phasing (LRP) or PLEP Angle Transfer Orbit 


terms of Kozai constant. TOF Lunes Fiyty to CLEP. 


Line of Nodes 
@ Lunar Flyby 


@ TESS PAM Burn 


Change since 
PDR Peer Review 
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cosAOP = ta - ay=e), 
er r 


Consistent with EXP-TESS- 
GSFC-RQMT-0015 Rev (-) 493 


L3 Mission Design Requirements 


ID Parent ID Title Requirement Compliance 

L3_ FD 21 MRD_54 | Launch Period FD shall design for at least 5 launch days in Comply. At least 
any given Lunar cycle. 9 sol’ns/mo for 

current period. 

L3_ FD_22 MRD_55 | Launch Window FD shall design for launch windows of at least Comply. Current 
5 minutes during each day of the launch period. | strategy meets 

req. 

L3_FD_27 MRD_42 | Commissioning FD shall design the phasing loops and post Comply. PAM at 

Duration lunar encounter transfer orbit to achieve < 43 days. 
mission orbit within 2 months after launch. 

L3_ FD_24 MRD_85 | Sun in Instrument FD shall design the PAM to occur when the sun | Comply. Basis for 

Boresight is not within a FOV of 54°x126° centered on sol’n selection. 
the camera boresight axis (X-Z plane) for 215 
minutes. 

L3_FD_28 | MRD_104 | Delta-V Budget FD shall design ascent-to-mission orbit to Comply. See 
require no more than 215 m/s delta-V with 99% | detailed analysis. 
probability of success. 

L3_FD_25 | MRD_129 | Maneuver Magnitude The largest maneuver magnitude shall be Comply. PAM < 
<95m/s. 75 m/s 

L3_FD 4 MRD_56 | Eclipse Frequency and | FD shall target a mission sequence that limits Comply. Ne-mere 

Duration the total number of eclipses from LV separation | #aAa++<4hr++ 
through the end of the prime mission to 2 <€hr 
eclipses with a maximum eclipse duration of 5 Needs updating 
irements added to flow frorn L2 hours, and 14 additional eclipses with a 
maximum eclipse duration of 4 hours. 


an 


ge since PDR Peer 
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eview 
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Lunar Flyby Orbit Geometry Options 


e Data shows 
generally best 


results for: iesconding icrandinn 
= Pre-flyby inbound outbound” ? inbound” 
= Post-flyby 
descending 
= Post-flyby Ascending node at encounter 
Z 
outbound 
e Pre-flyby 
ascending/ 
descending can “Descending “Descending 
be selected outbound” inbound” 


Figure 6: The four possible paths of the spacecraft following the lunar flyby. 


e For operational 
simplicity, we 
currently use 
ascending case 
only. 


e Implies short- 
_ coast solution at 
mEarth departure 


From Gangestad, J. et al. “A High Lunar Resonant Orbit for Lower Cost 
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Acronyms 


e CSALT Collocation Stand Alone Library and Toolkit 
e EMTG Evolutionary Mission Trajectory Generator 
e GMAT General Mission Analysis Tool 

e OGA Other Governmental Organization 

e PAM Post Apogee Maneuver 

e PLEP Post Lunar Encounter Periapsis 

e SNOPT Sparse Nonlinear Optimizer 


e TESS Transiting Exoplanet Survey Satellite 
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